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Simple thermodynamics considers kinetic energy to be an extensive variable which is proportional 
to the number, N, of particles. We present a quantum state of N non-interacting particles for which 
the kinetic energy increases quadratically with N. This enhancement effect is tied to the quantum 
centrifugal potential whose strength is quadratic in the number of dimensions of configuration space. 
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CN ■ I. INTRODUCTION 

The intensity of light radiated from independent dipoles is the number, N, of dipoles times the intensity of a single 
one. However, when the dipoles are located within a wavelength of the radiation the total intensity is iV 2 times the 
intensity of a sing le dipole @]. This supcrradiance effect is due to constructive interference between the individual 
dipoles. In the present article, we propose an enhancement effect of a similar strength for the kinetic energy of matter 
waves confined to a region whose radial extension in hyperspace is essentially independent of the number of particles. 

According to simple thermodynamics, kinetic energy is an extensive variable, that is, in general the average kinetic 
energy of an ensemble of N particles is linear in N Q. A similar dependence holds for the kinetic energy of an 
1 ordinary Bose-Einstein condensate (BEC) of N particles. However, we now show that for a special quantum state of 
N non-relativistic particles the kinetic energy increases as N 2 . This quantum state is completely symmetric under 
■ exchange of the coordinates of the particles, but there is no interaction between the particles. 

For the enhancement effect to occur, all particles need to have the same mass, but they do not necessarily have to 
fSJ . be identical particles, Nevertheless, the state could also correspond to N identical particles. In this case, it would be 
' the state of N bosons which are strongly entangled. 

Dimensional enhancement of kinetic energies is due to the wave nature of the atoms. It results from the form of 
the Laplacian in D dimensions, giving rise to the quantum centrifugal potential. We illustrate this phenomenon using 
N non-relativistic particles of identical mass, M, in three space dimensions. Here, we concentrate on the motional 
degrees of freedom, but do not take into account the internal structure of the particles. Hence, we deal with a D = 3N 
dimensional configuration space, and the wave function W = \&(xi, X2, ■ ■ ■ ,%d) depends on D coordinates. 

The proposed effect is most conspicuous for s-states, that is, when the wave function depends on the hyperradius, 
^ r = (x\ J r x\ J r- ■ ■ + x 2 D ) 1 / 2 , only. In this case, the wave function is completely symmetric under exchange of coordinates 
of the particles, corresponding to a bosonic state. We shall only consider s-states in the present work. 

Our paper is organized as follows: In Sec. [Il| we lay the ground work for the calculation of the kinetic energy of 
an ensemble of particles, by introducing the concept of the radial wave function in hyperspace. We illustrate it by 



using two examples related to a BEC in an isotropic harmonic trap. In Sec. Ill, we then turn to the discussion of the 
operator of kinetic energy and cast it into a form which brings out the quantum centrifugal potential. The latter is 
proportional to the square of the number of dimensions. It is this potential that may cause the kinetic energy to be 



quadratic in the number of particles, as discussed in Section [TV. 

In Sec. 0, we evaluate the kinetic energy of three different radial wave functions. The two motivated by BEC, 
and denoted by uo and u\, show a linear dependence on the number of particles. However, for the third radial wave 
function U2, which is independent of dimensions, we find a kinetic energy that depends on the square of the number 



of particl es. W e dedicate Section VI to a discussion of the origin of this enhancement effect. 

In Sec. VII, we then turn to a discussion of the quantum dynamics, starting from the wave functions uq, u%, and 
U2- We evaluate the time dependence of the average radial momentum following from these initial conditions. We 
note that the momentum corresponding to the BEC wave functions increases with a steepness that is proportional to 
the square root of the number of particles. In contrast, for the wave function ui we obtai n the remarkable result, that 



the steepness depends quadratically on the number of particles. We conclude in Section VIII with a brief summary. 
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II. WAVE FUNCTIONS 



For an s-state, the motion of the N particles is described by the normalized wave function = ty(r) or, equivalently, 
by the radial wave function u(r) defined by the relation 



*(r) 



1 u(r) 

m — d~ 

dd r 2 



where Sd denotes the total solid angle in D dimensions, 



S D = 



2n D / 2 
T{D/2) 



and u(r) is normalized such that 



\u{r)\ 2 dr = 1. 



(1) 



(2) 



(3) 



We shall look look at three particular wave functions of the above type. The first one may be constructed from a 
normalized Gaussian wave function, 



(4) 



in one dimension, by forming the product function 



#o(r) = <p(xi)<p(x 2 ) ■ ■ ■ <p(x D ) = [ — ) e ^ r2 



(5) 



with r being the hyperradius. This wave function describes the state of a BEC of N = D/3 non-interacting particles 
in an isotropic magnetic trap at zero temperature Q . The value of n is determined by the harmonic potential of the 
trap. 

Our second TV-particle wave function, \&i(r), may be constructed by retaining the wave function tp(x) for D — 1 
space directions, while taking a wave function of the form x 2 ip{x) for the last direction. Symmetrization by forming 
the coherent sum X^»=i xf exp(—^K 2 r 2 ), and subsequent normalization gives in fact 



*i(r) 



2k (D+4)/2 



ir D / 4 y/D{D + 2) 
The radial wave functions associated with 'I'oM and ^i(r) are 



2 -I K 2 r 2 

-.r z e 2 K r . 



uo(r)^Af r ( - D - 1 ^ 2 e-i K ~ r 



± K 2 r 2 
2 K r 
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with the normalization factors Afg and A/i being 



(6) 



(7a) 
(7b) 



r(f) 



& 2 j M = 



r(f + 2) 



(8) 



We defer the presentation of the third iV-particle wave function to Sec. |VB| 



III. OPERATOR OF KINETIC ENERGY 

We now consider the kinetic energies corresponding to the two above wave functions. The kinetic-energy operator, 



T = ———A( d \ 
2M 



(9) 



2 



is determined by the Laplacian 0] 



+ £^1*. + ^, (10) 
dr 2 r dr r 2 ' 

in D dimensions. The operator A 2 , which involves derivatives with respect to the D — 1 angles of hyperspace, is 
proportional to the angular momentum operator in configuration space [gj. 
As discussed in Ref. || , we may also write 

v 2 h 2 A 2 

Here, we have introduced the radial momentum 

1 (T „ „ r\ hi d p-i , , 

Pr = 5 (- P + P- -) = --Drr-^-r 2 , 12 

2 Vr rJi r T— or 



and the quantum fictitious potential 



h 2 (£>-!)(£> -3) 

^ (r)= 2M 4^ • (13) 



We emphasize that the quantum fictitious potential also emerges directly from acting with the D-dimensional 
Laplacian of Eq. ( |l0|) on the wave- function \&(r) of Eq. ([I]). This yields, in fact, 

92 ' V-W-V <r) . (14) 



Qi* 2 AlT 2 

Since Vq is proportional to the square of Planck's constant, it is a quantum potential with no classical analogue. 
Moreover, Vq depends inversely on the square of the hyperradius. This feature reminds us of the classical centrifugal 
potential, which gives rise to the non-inertial centrifugal force. Indeed, for D > 4 the potential Vq given by Eq. jl^ ) 
is positive and thus corresponds to a repulsive force. This property suggests the name quantum centrifugal potential. 
In accordance with this, some authors |Q absorb it into the last term of Eq. (|ll|). This, however, conceals the fact 
that it is a genuine part of the radial kinetic energy. 

We note that for D = 2, the potential Vq is negative corresponding to an attractive force. This centripetal force is 
unique to two dimensions and counterintuitive to the classical notion of the centrifugal force always being repulsive. 
To capture this contradiction, we have coined the phrase quantum anti-centrifugal potential for the potential Vq in 
the case of D = 2. In Refs. || and || we have focused on consequences of this attractive potential. However, in the 
present work we concentrate on the repulsive case corresponding to D > 4. 

We recognize that one and three dimensions are also special: the potential Vq vanishes. 



IV. AVERAGE KINETIC ENERGY 



The average kinetic energy 

fi 2 r°° r 

(f) = / dr r D - x \ dtt D -TA^)* (15) 

2M Jo J 

of the particles described by \& involves the Laplacian, Eq. (|lo|), and integrations over the hyperradius r and the solid 
angle £Id in the ZJ-dimensional hyperspace. With ^ = \&(r), as in the present work, the angular part of the Laplacian, 
that is, the operator A 2 /r 2 does not contribute to the integral. Hence, the kinetic energy may be said to be purely 
radial in hyperspace, and therefore also in the state spaces of the individual particles. 

When we introduce the radial wave function u(r), the average kinetic energy defined by Eq. (jl^) takes the form 

(f)=T = T r + T v . (16) 

Here the contribution 



3 



Tr 



dr u*(r) 



2M dr 2 



u(r) 



(17) 



is the average value of the operator p 2 /2M. We shall refer to it as the para-radial kinetic energy. 
Furthermore, the contribution 



T v = I dr V Q (r)\u{r)Y 
lo 



(18) 



results from Vq. 



V. EXAMPLES 



The dimension D of configuration space enters Vq quadratically. Hence, the total kinetic energy resulting from TV 
could in principle be quadratic in D. We recall that in the case of N particles in three space dimensions, we deal 
with a D = 3N dimensional configuration space. Consequently, for TV 3> 1, the strength S = (D — 1)(D — 3) = 
(37V — l)(3iV — 3) ~ 9iV 2 of the quantum non-inertial potential, is quadratic in the number of particles. According to 
Eq. (|l8|), the same thing may hold for the contribution Ty to the kinetic energy. 



A. Thermodynamic case 



However, this feature strongly depends on the form and, in particular, on the D-dependence of the radial wave 
function u. For example, the wave functions uq and u\ given by Eqs. ([7a]) and (fn]), respectively, yield the radial 
kinetic energies 



T 



(0) 



1 



1 



2(13- 2) 

and the quantum fictitious potential energies 



1 



2(L> + 2) 



T (0) _ 

J-v — 


(f 


- 1 - 


1 1 ) 


T (l) _ 




-3 + 


15 1 
2 D + 2 



(19) 

(20a) 
(20b) 



Here, we have introduced the kinetic energy e = (Tin) 2 / (2M) . 
Hence, the total kinetic energies 



T (0) 



D 



D 



D 



(21) 



are linear in D. In addition, the expression for involves correction terms independent of and inversely proportional 
to D. 

In the limit of a large number of particles, that is D — 3N 1, we recover the thermodynamic result 



T ~±De=-Ne 
2 2 



(22) 



for both states. 

It is int erest ing to note that the linear dependence on D shown in Eq. (|2l]) is due to the contribution Ty, Eqs. 
( |20a| ) and (20b), which emerges from the quantum centrifugal potential. The contribution from the para-radial kinetic 
energy is, according to Eq. (19), only weakly dependent on D. 
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B. Enhancement case 



We now introduce the new wave function 

u 2 (r) = M 2 exp 

with the normalization constant 



i (a 

2 



(23) 



M 2 = 09k) -* 



2Kx(2^/Pk) 



Here. 



/•OO 




/ r n exp 




/o 



r 



(24) 



(25) 



denotes the modified Bessel function |l0| of order 



In contrast to uo(r) and u\(r), Eqs. (7a) and (|7b|) , the wave function U2(r) is independent of the dimension D, 
that is, independent of the number of particles As a consequence, the cancellation of one power of D in the 

contribution Ty to the kinetic energy which appears for uq and U\, cannot take place for u%. Indeed, when we 
substitute the wave function u 2 into the definitions Eqs. (|T7]) and ([l8]) of the energies T r and Ty, and perform the 
integrations, we arrive at 



T (2) 



and 



T 



(2) 



1 K 2 (24fh€) 



(£>-!)(£> -3) 

4/3k 



(26) 



(27) 



(2) i ' i 

The para-radial contribution Tf ' is independent of D. In contrast, Ty' involves D, and hence the number N of 
particles quadratically. This enhancement results from the quantum centrifugal potential and reflects the constraint, 
that as we squeeze more particles into the state we do not alter the radial wave function u 2 . Indeed, u 2 is independent 
of D, and thus independent of N. Forcing additional particles into this state leads to a strong increase in energy. 

This situation is, to some extent, analogous to the problem of confining charged particles to a given domain of 
space. Due to their Coulomb interaction, 



(2) 



V c (r) 



<7i<72 



(28) 



two charges q\ and q 2 of the same polarity repel each other. This repulsion is proportional to the product of the 
charges and inversely proportional to their separation. 

In the case of neutral particles, the quantum centrifugal potential of Eq. (|l^) is proportional to the product 
(D — 1)(D — 3) of the dimensions. Hence, the dimension of configuration space plays the role of the charge. However, 
the analogy between Vq and Vc breaks down in the dependence on the separation. Indeed, the hyperradius r enters 
the denominator in a square, whereas in the Coulomb potential the separation enters only linearly. 



VI. ORIGIN OF ENHANCEMENT 

(2) 

This discussion suggests that we can interpret the energy T v associated with the quantum non-incrtial potential 
as localization energy. The localization energy is defined as the energy necessary to localize a quantum particle, with 

(2) 

the original energy T r , in a domain of hyperspace whose radial extension is essentially independent of D. This 
energy obviously depends on the dimension of the space in which we want to localize the particle, because the volume 
of the unit sphere in D dimensions decreases strongly for large D. Indeed, the expression (||) for Sd> together with 
an integration over r, easily gives the following expression for the volume 
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V D 



DT(D/2) 



(29) 



of the unit sphere. For D > 5 this is a decreasing function of D. 

(2) 

It is interesting that in two dimensions T v becomes negative. Hence, we do not need to perform work in order to 
localize the particle in this case, but rather gain energy. 

It is instructive to compare this quadratic dependence of the kinetic energy on the number of particles to the same 
scaling property of the total energy of iV fermions in a harmonic trap |T^| of frequency il. Indeed, when we assume 
that the fermions occupy the first iV states of the trap the total energy reads 



3=0 V ' 



(30) 



This simple calculation shows that the JV 2 -dependence arises from the fact that the particles are trapped in an 
external classical potential with a linear energy spectrum. In the case of the N particles described by the wave 
function of Eq. (pq), there is no classical potential present. Here, the N 2 -effect results from the quantum centrifugal 
potential, and from the fact that the strength S of this potential depends quadratically on the number of dimensions 
of configuration space. 

Nevertheless, there exists a potential V2 = V 2 (r) for which the wave function u 2 , Eq. (|2^), is an energy eigenstate. 
Indeed, when we differentiate u 2 twice with respect to r we find that U2 satisfies the time independent Schrodinger 
equation 



d 2 u{r) 2M 



dr 2 



+ ^r\ E -V2{r)]u{r)=Q 



with energy E — and the potential 



V 2 (r) = 

K ' 2M 



1/3 2 [3k (3 /K-y 
4^4 ~ ^2 ~ p + {2) 



(31) 



(32) 



Since U2 is independent of D, also the potential V2 has to be independent of D. 

We recall from Eq. (|l7|) that d 2 u(r)/dr 2 determines the para-radial kinetic energy rather than the total kinetic 
energy. To get the Schrodinger equation for the wave function \E f 2 (^), as related to u 2 {r) through Eq. (|lf), we draw on 
the expressions (|TT]) and (111), and get: 



- _A< D >* a (r) + [V 2 (r) - Vq (r)] * a (r) = B* 2 (r). (33) 
The confining potential in hyperspace is accordingly V 2 (r) — Vg(r). 

VII. QUANTUM DYNAMICS 

So far, we have focused on quantum kinematics. Now we turn to a discussion of the dynamics resulting from the 
three radial wave function uo(r), u±, and u 2 , considered as initial wave functions when we switch off the confining 
potentials at time t = 0, and hence allow for free-particle motion of the systems. Here we concentrate on the case of 
many particles, that is D = 3N ^> 3. 

The dynamics of this ensemble of free particles is governed by the Schrodinger equation 



dt v ' ' 2M 1 ' ; 



(34) 



By exploiting the relation ( |l4| ) for "J we obtain the following expression for the time evolution of the radial wave 
function: 



1 ~dt U V>> = 



-2Md^ + VQ{r) 



u(r, t). 



(35) 



For short times, we may obtain a first approximation to u(r, t) by replacing Eq. (pq) with the equation 



G 



iK-u(r,t) = [W(r) + V Q (r)]u(r,0), 



(36) 



where 



W(r) 



u(r,0) 



' 2M dr 2 



u(r,0). 



(37) 



This quantity is similar to Bohm's quantum potential (l3J . Integration gives 



u(r, t) ~ exp 



[W(r)+V Q (r)]t\u(r,0). 



(38) 



The formula ( pq ) allows us to derive an analytical expression for the short-time behavior of the average radial 
momentum 



(Pr>(*) 



dr r 



D-l 



dn D ^*(r, t) Pr ^{r,t). 



(39) 



By means of the definition Eq. ( [l2|) of the radial momentum p r and the ansatz Eq. (Q) for the radial wave function, 
this expression reduces to 



(Pr>(«) 



dru*{r, t) 



i dr 



u(r, t). 



(40) 



When we substitute the approximate solution Eq. (38) into this formula, we observe that W(r) does not contribute 
to the integral because u(r, 0) is real valued. Hence, we arrive at 



(Pr)(*) 



drF Q (r)\u(r,0)\ 2 t, 



(41) 



where 



F Q {r) 



dVq 
dr 



h 2 (D - 1)(D — 3) 



2M 



2r 3 



(42) 



This quantity, determined by the derivative of the potential Vq, is the quantum centrifugal force. The fact that the 
expression for (p r )(t) is independent of W(r) implies that we obtain the correct result for the short-time behavior of 
(p r )(t) by neglecting the para-radial kinetic energy operator compared to the potential Vq in Eq. (|35|). The result 
of doing so is known as the Raman-Nath approximation |l4|] . We also note that the integral preceding t in Eq. ( |4l| ) 
is nothing but d(p r )(t)/dt evaluated at t = 0. The expression ( [fl| ) may therefore also be obtained by differentiating 
both sides of Eq. (40|) with respect to t, while applying Ehrenfest's theorem pl|. 

To sum up, we have found that the initial average radial momentum increases linearly in time. The slope of the 
increase is determined by the average centrifugal force of the initial state. 

Let us now consider the specific wave functions uq, u\. and u^. We start our discussion with the wave functions 
i*o and u\ of Eqs. ((7a|) and (jFE]), respectively. Substitution into the expression (^|) for the average radial momentum 
yields 



r(V 



— TiK, 



(Pr) (1 Ht) 



r(#) 

1 rf^+i) fet 



With the help of the asymptotic formula [Q 

T(az + b) 

we can evaluate the ratios 



2Tre- az (az) az+b - 1 / 2 , 



Tin. 



r(V) 
r(#) 



D 



1/2 



r(^) 



D 



3/2 



(43) 



(44) 



(45) 
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in the limit D 3> 3. For both wave functions, this yields 



(Pr)(t)~ V2D^jjh K . (46) 

Hence, in the case of the initial wave function u and U\, the slope of the momentum increase is governed by the 
square root of the number of dimensions p6[ . 

Next, we turn to the wave function 112, Eq. (^3|). In this case, the average radial momentum Eq. (^) takes the form 

, ... (D - 1)(D - 3) K 2 (2y/PK) { et\ + 



Thus, the momentum corresponding to the wave function ui increases as D 2 , and hence depends quadratically on the 
number of particles. 

The explosion of the particles may be compared with the phenomenon of a Coulomb explosion . Atoms passing 
through a foil get stripped of some their electrons and become positively charged ions. Due to their Coulomb 
interaction, they repel each other and fly apart. In our case of neutral particles, it is the quantum centrifugal 
potential Vq, Eq. (O), which causes the explosion. 



VIII. SUMMARY 



We conclude by summarizing our main results. We have analyzed the kinetic energy and the dynamics of N non- 
interacting particles in free space. We have found a quantum state for which the average kinetic energy increases 
quadratically with the number of particles. Moreover, this kinetic energy stored in the state gets transferred into 
outgoing radial momentum in the time evaluation pursuant to the preparation of the wave function. In case of U2 
the explosion is more violent than for uq and u\, since the increase of momentum is proportional to N 2 rather than 
TV 1 / 2 . Both effects — the dimensional enhancement of kinetic energies and the iV 2 -explosion — are consequences of the 
quantum centrifugal potential. 
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